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^ ' Abstract 

^-H ■ This paper is concerned with the Schrodinger equation for atoms 

^~~^ \ and ions with A^ = 1 to 10 electrons. In the asymptotic hmit of large 

P_^' nuclear charge Z, we determine explicitly the low-lying energy levels 

i^ . and eigenstates. 

II I The asymptotic energies and wavefunctions are in good quantita- 

'^ ' tive agreement with experimental data for positive ions, and in ex- 

C^ . cellent qualitative agreement even for neutral atoms {Z = N). In 

G ! particular, the predicted ground state spin and angular momentum 

quantum numbers {^S for He, Be, Ne, ^S for H and Li, '^S for N, ^P 

CN ■ for B and F, and ^P for C and O) agree with experiment in every case. 

^ . The asymptotic Schrodinger ground states agree, up to small cor- 

^\i I rections, with the semi-empirical hydrogen orbital configurations de- 

^O ■ veloped by Bohr, Hund and Slater to explain the periodic table. In 

^-% . rare cases where our results deviate from this picture, such as the or- 

^^ ' dering of the lowest ^D" and '^S° states of the Carbon isoelectronic 

QQ • sequence, experiment confirms our, not Hund's, predictions. 

o: 

•^ ■ 1 Introduction 

c^ I How do the striking chemical differences between some elements, and the 

similarities between others, emerge from the universal laws of quantum me- 
chanics? In the physics and chemistry literature, this fundamental question is 
discussed on a semi-empirical level, via the "hydrogen orbital configurations" 
develo ped by Bohr, H und and Slater (see e.g. |Boh22t IHun251 [Sk3n| [CS351 
ILL77t ISchOll lAdPOlj ). or via numerical simulation of simplified quantum 
mechanical models (see e.g. |Har28t IHar57t iFFTTl ITTST94[ IJohnSj [BT86] ). 
In this article, we address this question from a mathematical perspective. 

The theoretical possibility of making chemically specific predictions was 
realized almost immediately after the Schrodinger equation had been intro- 
duced (see e.g. |Har28[ IDir29] ). But we are not aware of previous rigorous 



results, as mathematical research on the basic quantum mechanical equations 
has hitherto focused overwhelmingly on universal, qualitative properties. 

The perhaps most basic non-universal properties of atoms relevant to 
chemical behaviour are the total spin and angular momentum quantum num- 
bers S and L of the ground state, which describe the amount of symmetry 
under spin and spatial rotation. These two numbers not only determine the 
ground state dimension d, but, as argued below, they also allow to predict, 
up to at most two possibihties, the group of the atom in the periodic table. 

Other quantities of interest include the energy levels En and, more impor- 
tantly, energy differences such as spectral gaps En — Em (which govern the 
photon frequencies which the atom can emit or absorb) and the ionization 
energy (whose striking empirical periodicities lay at the origin of the design 
of the periodic table). 

Our principal result is that such quantitites can be extracted analytically 
from the many-electron Schrodinger equation in a natural scaling limit. More 
precisely, we show that for ions with A^ = 1 to 10 electrons, as the nuclear 
charge Z gets large the low-lying energy levels and eigenstates converge to 
well defined limits, which can be determined explicitly. In particular, this 
yields rigorous values of L, S, d for the ground state for all sufficiently large 
Z. See Theorems EZl El O and Tables [H El [ISl [H We call the above 
fixed- A^, large- Z limit is o- electronic limit, because it is realized physically by 
an iso-electronic sequence such as Li, Be"^, B"*""*", ... Note that this limit is 
different from the Thomas- Fermi limit N = Z ^ oo, which is of interest in 
other contexts but does not retain any chemical specificity. 

The asymptotic levels and eigenstates are in good quantitative agreement 
with experimental data for positive ions, and in excellent qualitative agree- 
ment even for neutral atoms {Z = N). In particular the predicted values of 
L, S and d (see Tabled]) agree with the experimental atomic values |RJK"'"07j 
in all cases. 
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Table 1: Angular momentum and spin quantum numbers and dimension of 
the Schrodinger ground state for large Z, as calculated in this paper. All 
numbers agree with the experimental values for neutral atoms |RJK"'"07] . 



The asymptotic ground states we calculate (see Theorem 13.11) provide for 



the first time a mathematical justification of the celebrated semi-empirical 
"hydrogen orbital configurations" developed notably by Bohr, Hund and 
Slater to explain the periodic table. In our approach, none of the underlying 
nontrivial postulates (electrons filling hydrogen orbitals, shell and sub-shell 
formation, sub-shell ordering rules such as 2s<2p, Hund's rules) need to be 
invoked, but are seen to emerge in a natural way. The only corrections are 
as follows (see Sections l3.ll and [HI for a detailed discussion): 

(1) Alongside each Slater determinant built from admissible hydrogen 
orbitals, the ground state must contain its orbit under the symmetry group 
5*0(3) X SU{2) X Z2 of the many-electron Schrodinger equation. 

(2) For the three elements Be, B, C, a ten to twenty percent admixture 
of a particular "higher sub-shell" configuration is also present, an effect we 
term 2s^-2p^ resonance. 

(3) In rare cases, such as that of the lowest ^D" and ^S° states of Carbon, 
the ordering of excited states disagrees with the semi-empirical Hund's rules 
(with experiment confirming our orderings). 

We now outline our mathematical strategy to obtain explicit asymptotic 
energy levels and eigenstates, focusing for simplicity on the ground state. 

The first step is the derivation of a simplified model governing the asymp- 
totics. This can be done via a scaling argument plus standard perturbation 
theory, as follows. For fixed A^ and large Z, attraction of an electron by the 
nucleus dominates its interaction with the other electrons, so one expects the 
true ground state to be close to the ground state of the corresponding system 
with electron interaction turned off (which is known explicitly via hydrogen 
atom theory). After a little more thought, one realizes that this cannot be 
quite correct. The non-interacting ground state eigenspaces of atoms happen 
to be highly degenerate (see Table [3]), but the underlying symmetry is broken 
by the interaction, so the true ground state eigenspaces should converge only 
to particular subspaces of them. (Experimentally, this phenomenon is well 
known, from observed energy splittings.) Mathematically, we prove that the 
difference between the Schrodinger ground states and the ground states of the 
problem PHP^> = E^^, where P is the projector onto the non-interacting 
ground state but H is the full Hamiltonian (eq. (^ below), tend to zero. 
We call this simplified problem PT model, because it corresponds to (i) re- 
scaling the problem so as to make the ground state of the reduced problem 
Z-independent, (ii) applying degenerate first order perturbation theory, (iii) 
undoing the re-scaling. Physically it corresponds to resolving, within the 
non-interacting ground state eigenspace, the fully interacting problem. 

The second step is to determine the lowest PT eigenspace. This requires 
a careful analysis of the interplay between hydrogen orbital formation (pro- 
moted by the Laplacian and electron- nuclei interaction), antisymmetry, spin, 
and electron interaction. More technically, the following difficulties arise, 
(i) The non-interacting ground state, i.e. the state space of the PT model, 
is of somewhat daunting looking dimension, e.g. 70 in case of Carbon (see 



Table [3D. 

(ii) The PT Hamiltonian PHP is easy to write down abstractly (as we have 
just done), but unknown; one needs to devise a method to determine it ex- 
plicitly. 

(iii) The PT model is a strongly interacting many-body model. 
Difficulties (i) and (iii) are overcome via careful use of the symmetry group 
of the original equation and its representation theory in terms of many-body 
spin and angular momentum operators, which allows one to split the Hamilto- 
nian PHP into small invariant blocks, (ii) is addressed by combining ideas 
from quantum chemistry which have not hitherto played a role in mathe- 
matical studies, such as Slater's rules |S096j (which allow to express the 
components of the Hamiltonian via six-dimensional integrals of a product of 
four hydrogen eigenstates and a Coulomb repulsion term), Fourier analysis 
(while the Fourier transform of individual hydrogen eigenstates is well known, 
here one requires the Fourier transform of pointwise products of these), and 
residue calculus. In principle, our methods apply to arbitrary atoms, except 
that the relevant PHP matrices can become significantly higher dimensional. 

One curious mathematical phenomenon we observe is that the Hamilto- 
nian PHP arising in the Z — > oo limit of quantum mechanics is always a 
rational matrix, despite H being a somewhat complicated partial differential 
operator and P a "transcendental" projector (onto tensor products of scaled 
hydrogen eigenfunctions such as 7r~^/^e~'^'). 

The structure of this paper is as follows. In Section 2 we describe the 
basic quantum mechanical equations and their symmetry group. In Section 
3 we state the asymptotic limit of the Schrodinger ground states for Li to Ne 
(see Theorem 13. ip . In Sections 4 and 5, we justify the reduction to the PT 
model and determine explicitly its state space. Sections 6-7 contain more 
technical material: the explicit determination of the PT Hamiltonian and the 
derivation of Theorem 13. ![ as well as of the excited states and levels of the 
PT model. Finally, sections 13.11 and [H] compare our results to experimental 
data and to methods in the physics and chemistry literature. 

2 Schrodinger equation and mathematical def- 
inition of basic quantitites of chemical physics 

The exact (nonrelativistic, Born-Oppenheimer) time-independent Schrodinger 
equation for atoms and ions is 

Hm = E^, (1) 

where, for nuclear charge Z > and A^ electrons and in atomic units, 

N ^ 

1=1 ' *' l<i<j<iV ' J^ 



E e M, and 

^eLldR'xZ^f). (3) 

Here and below the Xi G M^ are electronic coordinates, Sj G Z2 = {±|} are 
spin coordinates, and L^ is the usual Hilbert space of A^-electron functions 
^ : (M^ X Z2)^ — > C which are square- integrable, 

J] |^(a;i,si,...,a;7v,S7v)r = ll^lT < 00, (4) 

and satisfy the antisymmetry principle that, for all i and j, 

Sf (^ . . . , X j , S j , . . . , Xj y Sj y . . . J Sf (^ . . . , Xj , Sj , . . . , X j , S j , . . . J . 1^ J 

Mathematically, ilf is a bounded below, self-adjoint operator with domain 
L^ n H'^, where if ^ is the usual Sobolev space of L^ functions with second 
weak derivatives belonging to L^ [KatSlj . 

We are interested in the mathematical derivation of a number of quanti- 
ties of basic physical and chemical interest, and begin by recalling how these 
are defined in terms of the Schrodinger equation ([T]). 

Definitions, 1 An energy level of an atom or ion is an eigenvalue of the 
corresponding operator H. An eigenstate of the atom or ion is an eigen- 
state of H (i.e. a nonzero solution \E' to ([T]) belonging to the domain of 
H). By Zhislin's theorem ( |Zhi60j . see [Fri03] for a short proof), for atoms 
(A^ = Z) and positive ions (A^ < Z) there exist countably many energy 
levels El < E2 < ... below the bottom of the essential spectrum of H, the 
corresponding eigenspaces being finite-dimensional. Ei is called the ground 
state energy and the corresponding eigenspace is known as the ground state. 
Eigenspaces corresponding to the higher energy levels are known as excited 
states. The excitation energy or spectral gap of an excited state with energy 
Em is defined to be Em~Ei. Physically it corresponds to the energy required 
to promote the electrons from the ground state to the excited state. 

Besides the quantized energy levels En, there exist important additional dis- 
crete quantum numbers associated with the atomic Schrodinger equation 
which arise from its symmetries. Their precise definition, albeit very natu- 
ral, takes a little more work. 

The model ([T]), ([2]), dl]), © is invariant under 

(i) simultaneous rotation of all electron positions about the origin, 
\E'(xi, si, .., xn, sn) >— > '^{R^xi, si, .., R^xn, sn), R G S0{3) 

(ii) simultaneous rotation of all electron spins 
(by a unitary matrix U G SU{2)) 

(iii) simultaneous inversion of all electron positions at the origin, 
^(xi, si, .., xn, Sn) ^ ^(-a;i, si, .., -xn, sn) =■ -R^- 



(In group theory language, the symmetry group is 50(3) x SU{2) x Z2, the 
third factor being the inversion group consisting of R and the identity. When 
N = 1, there exists an additional symmetry, which gives rise to conservation 
of the quantized Runge-Lenz vector; but it is broken by the interaction term 
in ([2]) when N > 2. Note also that even though the Hamiltonian ([21) is 
invariant under the larger group of non-simultaneous rotation of spins, the 
antisymmetry condition (^ is not.) 

The conserved quantities, i.e., operators which commute with the Hamil- 
tonian, which arise from the above symmetries are 

(i) L = Yli=iiL{j) (many-electron angular momentum operator) 

(ii) S = ^,=1 S(j) (many-electron spin operator) 

(iii) R (parity operator), 
where 

Uj) = I L,{j) I , S(j) 





and La(j), Sa(j) (a = 1,2,3) denote the usual angular momentum re- 
spectively spin operators acting on the j^'^ coordinate. Explicitly, on A^- 
electron states \E'(xi, si, .., xat, sat), Xj e M^, Sj G {±|}, and denoting Xj = 
(^y(i)^y(2)^y(3)-)^ Lq,(j) is the partial differential operator 

and Sa{j) is multiplication by a Pauli matrix, 

(S'„(j)^)(Xi,Si,...,Xj, i,...,XAr, Stv) \ _ / ^(Xi, Si, . . . , Xj, i, . . . , Z7V, S^v 



{Sa{j)'^){xi, Si, . . . , Xj, -^, . . . , Xn , Sn) J "V "^{Xl, Si, . . . , Xj, -I, . . . , Xn, Sn) 

where the a^ are the Pauli matrices 



'a 



1 /O 1\ 1 /O -A 1 fl 



^^ ^" 2 *^1 Oj' ''''■~2\t yl ' ''' - 2 \0 -1 

The fact that the operators (i), (ii), (iii) commute with the Hamiltonian 
([2]) can be checked by direct inspection using the above formulae, without 
reference to the underlying symmetry group. 

The components of total angular momentum and total spin. La = X]j=i -^a(j) 
and Sa = J2j=i ^a{j), obey the usual commutator relations 

[La, Lp] = iL^, [Sa, Sp] = iS^ (a, /5, 7 cychc). 

Angular momentum representation theory, together with simple considera- 
tions concerning the above specific action of the operators on iV-electron 
states, yields the following well known facts (see e.g. |FriXX] ) . 
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Lemma 2.1. (a) For arbitrary N and Z, a set of operators which commutes 
with the Hamiltonian H and with each other is given by 

L , L3, S , 5*3, R. (7) 

(b) The eigenvalues o/L^, S^, and R (acting on L^((M^ x ^2)^)) are, re- 
spectively, 



L{L + 1), L = 0,1,2,..., 






(8) 


s(s + i), s={|-i;|-; 


N 
••' 2 ' 

N 
•' 2 ' 


A^ odd, 
N even 


(9) 


p=±l. 






(10) 



(c) For fixed L, S and p, on any joint eigenspace of H , L^, S^ and R, 
L3 has eigenvalues M = —L,—L + 1,...,L, and S3 has eigenvalues Ms = 
—S, —S + 1, ..., S. In particular, the eigenspace has dimension greater or equal 
to (2L + 1) ■ {2S + 1), with equality in the case when the joint eigenspaces of 
H and the operators ^ are non- degenerate (i.e., one- dimensional) . 

Here we have employed the usual notation L^ = L^+Lg+Lg (and analogously 
forS'). 

From the above we see that the main differences between the symmetries 
of many-electron atoms and those of hydrogen are the absence of an analo- 
gon of the quantized Runge-Lenz vector and the nontrivial action of the spin 
operator S^ (for A^ = 1 it is equal to the trivial operator |/). 

Definitions, 2. The values of L, S and p for eigenstates of L^, S^, R are 
called the total angular momentum quantum number, the total spin quantum 
number, and the parity of the state. 

From the above lemma we see that for each energy level of H there ex- 
ist unique quantum numbers L, S, p characterizing the symmetry of the 
eigenspace (except in "non-generic" cases where the joint eigenspaces of H 
and the operators ([7]) are degenerate, in which case there exists a unique finite 
set of such quantum numbers). In the chemistry literature these numbers for 
an energy level are usually given in the form 25+1 j^v^ where L corresponds 
to X via — s> S", 1 — s> P, 2 — > D, 3 ^ -F, and where no superscript v means 
p = 1, and V = (for odd) stands for p = —1. For example, the Carbon 
values L = l, S = l, p=l from Table [H would be denoted ^P, and the 
Nitrogen values L = 0, S = 3/2, p = — 1 by ^5'°. 

Of particular physical and chemical interest is the energetic ordering in 
which different combinations of L and 5* appear in the spectrum of the Hamil- 
tonian ([2]); see Section [71 



3 Ground states of the first ten atoms in the 
large Z limit 

We are now in a position to state a principal result of this paper. 

Notation \ri\...ri^) denotes the Slater determinant (or antisymmetrized 
tensor product) of the orbitals ?7jGL^(R'^ x Z2), \ri\...ri^){x\^s\^ ..^Xfq^sj^^ = 
(N\)~^^'^ det{r]i(xj, Sj)f^j^i). t, | stands for the spin-up and spin-down 
orbitals (f){x)6±i/2{s) G L^(]R^ x Z2). For a linear operator on the iV-electron 
Hilbert space (^, \\\A\\\ denotes the usual operator norm sup{||y4\l/|| : \1/ G 
L2((R3xZ2n, ||vl/|| = l}. 

Theorem 3.1. Consider the atomic Schrodinger equation (QP, ^, ^ for 
the atom/ion with N =1 to 10 electrons and nuclear charge Z . 
(i) For sufficiently large Z , the ground state has the spin and angular mo- 
mentum quantum number S , L and the dimension given in TableUl 
(a) In the limit Z -^ 00, the ground state is asymptotic to the explicit vec- 
tor space given in Tahle\E, in the sense that the projection operators Pq, Pq 
onto these spaces satisfy lim^^oo lll-fo ~ -Polll = 0. Here 1, 2, 3, 4, 5 stand 
for the scaled hydrogen orbitals (mathematically: hypergeometric functions) 
0isT, 02sT, 02^3?, 02piT, 02p2T f^om ([^, (E^j, and T, 2, 3, 4, 5 denote the 
corresponding spin-down orbitals. 

In fact, all low-lying energy levels and eigenstates can be determined exactly 
in the above limit. See Theorem 17.11 Note also that, as we will see below, 
the asymptotic ground states in the table are the exact ground states of the 
limit model ( TT6l) . ( TT71) . and become Z-independent after the re-scaling ( TTTl) . 
The derivation of these results requires five steps: 

1. Reduction to the finite-dimensional problem PHP"^ = E^^, \E' G Vo(A^) 
described in (ITBD. (IT71). 

2. Explicit determination of Vo(A^). 

3. Choice of a suitable basis of Vo(A^) making PHP particularly simple. 

4. Explicit determination of the dxd (i.e., in case of Carbon, 70 x 70) matrix 
representing the Hamiltonian PHP from eq. ( fT6l) in this basis. 

5. Spectral analysis of this matrix. 

These steps are carried out in Sections HI [5], 16.31 I6.4H6.61 and [3 Steps 1. 
and 2. follow from standard perturbation theory respectively hydrogen atom 
theory. 3. is achieved by a basis adapted to the symmetries of PHP (see 
Lemma 12.11) leading to block diagonal structure. 4. exploits, in addition, 
the fact that the Hamiltonian contains only one-body and two-body terms, 
allowing to reduce evaluation of the required iV-electron matrix elements 
(^|iJ|\l/), which are integrals over M'^^, to 1- and 2-electron matrix elements. 



Iso-electronic 
sequence 



Symmetry 



Ground state 



Dimension 



"S 



|1): IT) 



He 



111) 



Li 



1 112), |1T2) 



Be 



|1122) + c^(|1133) + |1144) + |1155))) 
■ = -^(2^1509308377 - 69821) = -0.2310995 . 



\1122t}+cj^{\lUjj} + \Uikk} 

|lT22l)+c^(|lHj7) + |imfc) 
(i,j,fe) = (3,4,5),(4,5,3),(5,3,4) 
"33^(^733174301809 - 809747) = -0.1670823 . 



3p 



1 
/T+c 



^{\n22tj} + c\UkMj)) 

^(|1T2W) + |lT22*j)) + cj.^{\Ukk{J} + llU-Ej) 
1 



,(\1122tj} + c\llkkij}) 
(?:,j,A:) = (3,4,5), (4,5,3), (5,3,4) 
r(^221876564389 - 460642) = -0.1056317 . . . 



11122345) 
^(|1T22345) + |1T22345) + |lT22345)) 
^(|1T22345) + |1T22345) + |lT22345)) 

|1T22345) 



O 



|1122MJfc) 
^ ( I lT22il jk) + I lT22iIj fe) ) 

\1122iijk} 
(?:,j,A:) = (3,4,5),(4,5,3),(5,3,4) 



2po 



\n22fijjk} 

\n22iijjk) 
(i,j,fc) = (3,4,5),(4,5,3),(5,3,4) 



Ne 



11122334455) 



Table 2: Ground states of the atomic Schrodinger equation in the hmit Z -^ 
oo. The indicated wave functions are an ortho normal basis of the ground 
state. See Theorem 13.11 for notation. The symmetry agrees with experiment 
for each sequence and all Z. 



3.1 Comparison with the semi-empirical Bohr-Hund- 
Slater picture of the periodic table 

The result of Theorem 13. II provides a mathematical justification of the semi- 
empirical "Aufbau principle" (from the German word for building up) de- 
veloped notably by Bohr, Hund, and Slater to explain the periodic table 
|Boh221 IHun25t [LTTTTI ISchnit lAdPOlj . The aufbau principle is based on 
three semi-empirical postulates: 

(a) Each electron in an atom occupies a hydrogenic orbitaljj 



^In fact, in Bohr's and Hund's original works |Boh221 rHun25| . which narrowly predate 
the Schrodinger equation, the electrons were supposed to occupy hydrogenic Bohr orbits. 



(b) Sub-shell ordering The orbitals in each hydrogen energy level, or shell, 
form sub-shells which are occupied in the order Is 2s 2p 3s 3p 4s 3d . 

(c) Hund's rule Within any partially filled sub-shell, the electrons adopt a 
configuration with the greatest possible number of aligned spins. 

Thus, for example, in Carbon the six electrons would occupy the orbitals 
Is t Is I 2s t 2s I 2pi I 2p2 T (note that the alternative choices 2pi [ or 2p2 i 
for the last orbital would be consistent with (b) but not (c)). 

This beautiful heuristic picture is seen to emerge in Theorem 13.11 in a 
natural way, without reliance on the above nontrivial postulates or numerical 
simulations (up to small but interesting corrections). 

(i) For seven out of ten elements (H, He, Li, N, O, F, Ne), the Aufbau 
principle configuration (when interpreted not as a collection of individual 
electronic states, but as a Slater-determinantal many-electron wavefunction) 
is an element of the asymptotic Schrodinger ground state. 

(ii) For the remaining three elements (Be, B, C), the Aufbau principle con- 
figuration is the dominant part of an element of the asymptotic Schrodinger 
ground state. 

The corrections to the semi-empirical rules emerging in the large Z limit 
of quantum mechanics are: 

(1) Alongside any Slater determinant, the asymptotic ground state contains 
its orbit under the symmetry group of the many-electron Schrodinger equa- 
tion (see Lemma [2. II) . 

(2) The corrections to the Aufbau principle configurations in (ii) come from 
different sub-shells, indicating that rule (b) is not strictly obeyed. We term 
this effect, which does not seem to have received mathematical attention, 
2s^— 2p^ resonances: besides the aufbau principle configurations ls^2s^, 
ls^2s^2p, and ls^2s^2p^, a significant percentage is also present of, respec- 
tively, ls^2p^, ls^2p^, and ls^2p^. This could be described as a resonance of 
the standard configuration with a configuration in which the two 2s orbitals 
have been substituted by two 2p orbitals. Why this phenonenon occurs only 
for Be, B, C has a simple group-theoretic reason: Tables [12], [13] show that 
such a substitution which preserves the total quantum numbers L and S is 
only possible in these three cases. 

(3) For excited states, Hund's rules are in rare cases found to disagree with 
the experimental and mathematical results; see Section [8] 

3.2 Comparison with other approaches 

First we comment on the asymptotic regime in which the above picture 
emerges, namely N fixed, Z -^ oo. These parameters are the only ones con- 
tained in the electronic Schrodinger equation that can be varied in ground 
states of a physical system. To capture chemical specificity, A^ must be kept 
fixed. A priori, Z could be made either small or large; but making it small 
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leads to non-existence of bound electrons [Lie84j . This leaves Z ^ cxd, i.e. 
the limit studied above, as the only option. 

This limit has attracted considerable previous attention in the physics 
and chemistry literature in connection with asymptotic expansions of energies 
and numerical evaluation of the coefficients. See Section 4 and the references 
given there. But we do not know of any effort to derive a numerical analogue 
of Table [2] in this way (which would correspond to truncating the expansion 
at first order in l/Z, see Theorem 14.11) . even though, in principle, the tools 
to do so were available. In fact, in the only case we are aware of where first 
order wavefunctions are given |Lay59[ Table 1, p. 288], those for Be, B, C 



are incorrect (the standard aufbau principle configurations are given, instead 
of the correlated states in Table [2]). 

Another interesting strategy would be to vary mathematical parameters 
in the Schrodinger equation which cannot be varied physically, such as h ^ 
(semiclassical limit), ot D ^ oo where D is the dimensionality of single- 
particle state space. Despite interesting results (see e.g. the semiclassical 
picture of highly excited states of Helium in [TRROOj , and the total atomic 
energies in |Loe86j via large-D asymptotics |Wit80t[Her86j ). these ideas have 
so far not been aimed at, or led to, explanations of the main features of the 
periodic table. 

The large majority of the literature on atomic systems is computational, 
and does not take the Schrodinger equation, but various simplified models 
as a starting point. Explanations of the filling order (postulate (b) above) in 
terms of numerical solutions of the Hartree and Hartree-Fock equations began 
with the pionieering work of Hartree on Rubidium |Har28] : for treatments 
of large classes of atoms see e.g. |Har57l IFF77t ITTST94] . Note that these 
models assume a significant part of postulate (a) from the outset, namely 
that electrons occupy individual orbitals and that these have Is, 2s, 2p, ... 
symmetry. The asymptotic Schrodinger ground states which we determined 
above beautifully illustrate both the power of the Slater determinant ansatz 
in Hartree-Fock theory and its limitations: in several cases other than noble 
gases the ground state contains a determinantal state, but for some of the 
atoms it does not. 

Finally we remark that the refined, multi-configurational methods de- 
veloped in quantum chemistry (see e.g. |FF77t IBT861 IS096] ) assume, on 
a somewhat ad hoc basis, exactly the structure of the wavefunctions which 
emerges naturally in Table [21 namely finite linear combinations of Slater 
determinants. The only reason preventing these methods from being asymp- 
totically exact within numerical error as Z — > oo is the use of Gaussians to 
represent the orbitals (see our companion paper |FG09] ). 
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4 Reduction to Perturbation Theory Model 

The first step in establishing the above result is to show that ([T]) simphfies to 
a finite-dimensional model in the limit of fixed electron number A^ and large 
nuclear charge Z. 

If \l/ solves the original Schrodinger equation ([T]), then an elementary 
calculation shows that its rescaling 



^(Xi, Si, ... , X,v, STv) = Z-=^^/2^(Z-Vi, Si, . . 


. , Z XAr,SAr) 


(11) 


solves the equation 




(12) 


where Rq is the Z-independent Hamiltonian 






1=1 ' 




(13) 



(acting on the A^-electron Hilbert space ([3])) and E = -^E. The elementary 
but important observation now is that the interaction term -^Vee in (TT^ 
becomes small if Z is large, allowiong to treat the interaction by perturbation 
theory. 

Let us first derive the ensuing perturbation-theoretic model informally, 
then formulate a theorem. By first order Rayleigh-Schrodinger perturbation 
theory, flT^ is expected to be well approximated by 

P{Ho + ^Vee)P^ = E^, ^ eVo, P = orth.projector onto Vq (14) 
Vo = ground state eigenspace of Hq, (15) 

with Hq as in flT^ . Now we undo the rescaling fITT]) . This yields the model 

PHP-^ = E^, ^ G Fo, P = orth. projector onto Vq, (16) 

Vo = ground state eigenspace of Hq, Hq = N^ ( A^;- — -j — - )(17) 

i=l ' *' 

where H is the original Hamiltonian ([2]). 

We call eqns. (fT6l) . (fTTl) the PT model. While it is still a fully interacting 
quantum many-body model, the key simplification is that the space Vq is 
finite-dimensional. Its dimension for different atoms is easily read off from 
Lemma 15.11 below: 

An important feature of the PT model is that it retains the full symmetries 
of the atomic Schrodinger equation. 

Lemma 4.1. For arbitrary N and Z, with P as defined above and with H 

denoting the Hamiltonian ^, the operators ^ 

(i) leave the ground state Vq of Hq invariant 

(a) commute with the PT Hamiltonian PHP : Vq -^ Vq. 

12 



Atom 


He 


Li 


Be 


B 


C 


N 





F 


Ne 


dim Vo 


2 
1 


3 

8 


4 
28 


5 
56 


6 
70 


7 
56 


8 
28 


9 

8 


10 
1 



Table 3: Dimensions of degenerate Hq ground states, as given by Lemma EIH 



Proof By direct inspection the operators (171) commute with Hq. Since Vq is 
an eigenspace of Ho, they must therefore leave Vq invariant, and commute 
with the projector P onto Vq. As already shown (see Lemma [2. ip . the oper- 
ators ([7j) also commute with H, and hence with the composition PHP. 

We now come to the rigorous justification of the PT model (iT6l) . (ITTl) . 

Theorem 4.1. Let A^ = 1, . . . , 10, Z > 0, and let n{N) be the number of 
energy levels of the PT model ^E), p7| j. Then: 

(a) For all sufficiently large Z, the lowest n{N) energy levels Ei{N,Z) < 
■ ■ ■ < En{N) {N, Z) of the full Hamiltonian ^ have exactly the same dimen- 
sion, total spin quantum number, total angular momentum quantum num- 
ber, and parity as the corresponding PT energy levels E['^{N, Z) < ■ ■ ■ < 

(b) The lowest n{N) energy levels of the full Hamiltonian have the asymptotic 
expansion 

E,{N,Z) _Ef^{N,Z) ^,,, 1^, 



Z2 



Z2 



+ 






E(o) 



Z^^ 



(1) 



+ 



■ 1 ' 



as Z 



oo, 



where E^^^ is the lowest eigenvalue of Hq and the Ej are the energy levels 

ofPVeeP OnVo. 

(c) The projectors Pi, ... , Pn{N) onto the lowest n{N) eigenspaces of the full 
Hamiltonian satisfy 



yPT\ 



as Z —^ oo, 



where the P^^ are the corresponding projectors for the PT model. 

The idea that for large Z the inter-electron term E- provides the first order 
correction to the non-interacting energy is well known in the physics literature 
(see e.g. |Hyl30i IBSSTI 15^621 ISW671 IRD711 IWil84j who treat non-degenerate 
eigenvalues, and see |Lay59] who gives an expansion similar to that in (b), not 
accounting for multiplicities, and numerical tables of E, in the degenerate 
case) . The main new insight here is the absence of further splittings at higher 
orders of perturbation theory (see statement (a) in the theorem). This is 
remarkable, considering that it fails in the simple 3x3 matrix example 

^0 1 0\ /I 0^ 

Hit) =11 0+eO -1 
,0 1/ \0 0, 
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The eigenvalues are \/l + e^, 1, — Vl + e^, and hence nondegenerate for e 7^ 0, 
but the leading eigenvalue is degenerate in first order perturbation theory. 

Proof Let Ef"^ < ... < E^Tj^-, be the energy levels of the scaled model 

dH]), ( JTSj) . and let dj, Pf^ be the corresponding eigenspace dimensions and 
eigenspace projectors. By perturbation theory for relatively bounded pertur- 
bations of self-adjoint operators (see e.g. |Kat95l[FriXXj ). exactly dj eigenval- 
ues of the scaled Schrodinger equation (fT2ll including multiplicity are asymp- 
totic to first order in 1/Z to the j^^ eigenvalue of ( iT6l) . ( ITTll . More precisely: 
The lowest di + ... + dn{N) eigenvalues of (iT2l) including multiplicity, labelled 
E,-fc, J = l,..,n(Ar), k = l,..,ci„ Ei,i < ... < Ei^d, < E2,i < ... < E2,d, < ..., 
satisfy 

4fc = EW + i^f +0(^) asZ^oo, k=l,..,d^. 

Moreover the projector Pj onto the span of these dj eigenstates satisfies 

|||Pj-Pf^||| ^0 asZ^cx). (18) 

Next, we investigate the Schrodinger eigenspace dimensions. By Lemma HIT], 
each PT eigenspace possesses well defined spin, angular momentum and par- 
ity quantum numbers L, S, and p, and by inspection of the explicit formulae 
in Theorem 17.11 below, the space has minimal dimension subject to these 
numbers. On the other hand, by (ITSll . for sufficiently large Z these numbers 
must agree with those of the eigenspaces of (IT2l) : hence by Lemma [2. II (c), 
Ej^i = ... = Ej^d ■ Note that without the information on minimality of the 
PT dimensions, we would not be able to exclude the possibility of further 
splittings of the Schrodinger eigenvalues beyond the PT splittings, at higher 
orders of perturbation theory; this is the only reason why the restriction 
A^ < 10 is needed. 

The theorem now follows by applying the isometric scaling transformation 

5 State space of the PT model 

The important starting point for solving the PT model is the fact that its 
state space, the GS of Hq, can be determined explicitly. This will follow 
from the exact solubility of the Schrodinger equation of hydrogen and basic 
many-body arguments. To explain these matters, we start from the hydrogen 
atom Hamiltonian 

/i = -|A-^, (19) 

X G M.^, acting on L^(]R^ x Z2). For hydrogen, Z = 1, but the parameter 
Z > will be useful later. Its eigenvalues are given by (see e.g. |Gri95] ) 

e. = -£^, neN, (20) 
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and have corresponding 2n^-diniensional eigenspaces with orthonormal basis 

Bn = {(pnlmsix, s) \ £ = 0, . . . , U-l, TTl = -£, -i+1, ...,£, S = -^, ^} , (21) 

where (pnems ^ -L^(M'^ x Z2) is the, up to normahzation, unique eigenfunction 
of h, L^, L3 and ^3 with eigenvalues — Z^/(2n^), £{£+!), m and s respectively. 
For later it will be useful to have an explicit form for these so-called 
hydrogen orbitals, which in polar coordinates with spin coordinate s G Z2 
are given by 

0n,Ln^,.(r, 9, 0, s) = (P^.^^T, 9, 0)(5<,(s) = Z''/^ R^^eiZr)Ye,UO, <f>Kis) (22) 
where 

«-«-((!)^ll(;^)"^--(l)'^S'.(|). (23, 

Here L^{x) is a generalized Laguerre polynomial and Y£^m{9, (p) is a spherical 
harmonic |AS72] . In cartesian coordinates, the n = 1 and n = 2 orbitals are 

^3/2 

0i,o,o(a;) = ^^e~^l''l =: 0is(x), 



^2,0,0(a:) = 4^ f 1 - ^ ) e-^l^l/^ =: ^2s[X) 



V ^TT V 2 

2-5/2 1^1 

<^2,i,o(a;) = -j==X3e-^^ =: 02p3, (24) 

V 0277 

J, / \ ^^^^ ^1 =^^^2 -Zbl/2 i. 

Often, it is convenient to work - instead of the last two functions - with their 
real linear combinations 

—=a:,e-^-=:hpA^), J = 1,2. (25) 

The following lemma describes how the eigenfunctions for the non-interacting 
many-electron system are formed from these one-electron eigenfunctions. 

Lemma 5.1. (Standard "folklore", see JFriXX^ for a rigorous proof) (a) 
The lowest eigenvalue of the operator 

N N 

4 = 1 j = l ' ' 

on the space L^(^(R^ x Z2)^) of square-integrable functions \E' : (M^ x Z2)^ -^ 
C satisfying the antisymmetry condition ^ is E = Y^^^^en, where ei<e2< 
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... is an ordered list including multiplicity of the hydrogen eigenvalues l[2U]) . 
(b) The corresponding eigenspace is 

Vo = Span I Ixi . . . Xd^i^h ■ ■ ■ i'iN-^d, ) 

l<^l<■■■<^N-d,<2{n, + lY}y (26) 

where the functions Xi andipi ("core orbitals" and "valence orhitals") and the 
integers d* and n^ ("number of core orbitals" and "number of closed shells") 
are defined as follows: d^{N) is the largest number of form X]?=i 2j^ which 
is less or equal to N , n^ (N) is the corresponding value of n, 

(union of the ON bases l[21\) of the first n* hydrogen eigenspaces), and 

(ON basis of the (ra* + 1)** hydrogen eigenspace) . 

Thus the ground state of the non-interacting Haniiltonian is spanned by 
Slater determinants (ahas antisymmetrized tensor products) formed from 
scaled hydrogen orbitals, "filled" in order of increasing one-electron energy. 

Due to the freedom of choosing any N — d^, eigenfunctions ipt ( "valence 
orbitals" ) from the basis of the highest relevant hydrogen eigenspace, whose 
dimension is 2(n^, + 1)^, the noninteracting GS has typically a large degen- 
eracy: 



do:^d.mGSofi7o-^^^:— J. (27) 

Specialization to the second row atoms and their isoelectronic ions 

(N=3,...,10) In this case, the number d^ of core orbitals equals 2, the number 
n^, of closed shells equals 1, and the dimension 2{n^, + 1)^ of the hydrogen 
eigenspace from which the valence orbitals are selected equals 8. Thus by 
( 1271) . the dimension of the ground state equals 



dt, 

\N-2 

These numbers are given in Table El The set of core respectively valence 
orbitals is (using the real orbitals (p2pi, 4>2p2 instead of 02p±) 

{Xi,X2} = {0isT, 0i4}, (28) 

{^l,...,tps} = {02. T, <p2s i, </'2pi T, hpi L hp2 T, <f)2p2 i, </>2p3 T, hps I}- (29) 

Here we have employed the standard notation 1, 1 for the two spin orbitals 

0(x)(5±i/2((t). 

Finally, the ground state of Hq is 
Vo{N) = Span ||xiX2V'n • • • ^1^-2) 1 1 < «i < ■ ■ • < tN-2 < s}. (30) 
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6 Determining the matrix PHP 

In this section we determine explicitly the Hamiltonian matrices PHP, for 
all second period atoms. 

Most of our arguments do not rely on the special radial form of the hy- 
drogen orbitals fl^ . fl2I?l) appearing in the definition of the subspace Vo(A^). 
Hence in this section, unless stated otherwise, Vo(A^) denotes the space 
( !28l) . (!29l) with the more general orbitals 



<^is{x) = Ri{\x\), (p2s{x) = R2{\x\), V52p,(a;) = -RalkDxi (i = 1,2,3), (31) 

where the (/j's are in L^(M^) with norm one, i?j : M — > M, and J^ Ri{r)R2{r)r'^dr 
0. 



6.1 Spin and angular momentum calculus on Slater 
determinants 

The action of the spin and angular momentum operators on Vo(A^) can be 
calculated from their action on the orbitals fl2^ and fpU|) together with the 
following simple identitites for the action of linear operators of form 

N N 

i=l ij=l 

on Slater determinants, where 6 is a linear operator on //^(R'^ x Z2): 

N 

B Ixi, • • • , Xiv) = X] 1^1' • • • ' ^^^' • • • ' Xiv), (32) 

i=l 

N 

j=l 
+ 2 ^ |xi,...,&Xi,---,&Xi,---,X7v). (33) 

l<i<j<N 

Direct calculations show that, for any two spatial orbitals '?/', € L^(M^), 
and orthogonal spin states a, /3 : Z2 — > C, 

S_ ■ Sjpa = \ipa, 
(i5(l) ■ Si2))^lja ®(t)f5= ltp(3 ®(pa- \^a ® 0/3, 
{S{1) ■ S_{2))i)a ®(f)a = \%l)a ® (pa. 

In particular, by (!32|) and (!33l) . S^lipaipP) = and S^\ipo;^P) = (| + 



^4 

inr'^YM/ " 41 



3.)|^«V^/?)-2-f|V^«^/?) = 0. 
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The angular momentum operators ([6]) act on the orbitals fl3Tl) as follows, 
independently of the choice of spin a : Z2 ^ C: 

LjifiisO = Lj(p2sa = Ljip2p,a = 0, 
Lj+iip2p^a = - iip2p,_^a, Lj^np2p^a = iip2p,+ia, j = 1, 2, 3, 

where the indices are understood modulo three. Hence we need only consider 
the action of L^ on ^2pi, giving, for any two spin states a and P, and i y^ j, 

L- L(p2p,<y = 2v52p,a 
(L(l) • L{2))i^2p,a ® <^2p,f3 = -{^2n„^a ® <^2p,.,P + '^2p,+^a ® (^2p,+i/5) 
(L(l) ■ L(2))^2p,a ® ^2pj(3 = ^2p,a ® ip2p,f3. 

Finally we see that, for any spin state a, Ri^ua = ipua, Rip2sO: = (p2sO: 
and R(p2p,a = -v?2p,a, i = 1, 2, 3. 

A useful and well known consequence of the above is that the pair of Is 
orbitals makes no contribution to spin, angular momentum, and parity on 
the space ( l30l) . More precisely: 

Lemma 6.1. (See [FriOX]) The matrix of any of the operators ^ on Vo(A^) 
with respect to the basis (SD^) is the same as that on the corresponding fewer- 
particle space obtained by deleting the orbitals xi; X2, with respect to the 
corresponding basis {|^j^ . . . V'i]v-2) 1 1 < "^i < ■ ■ ■ < in -2 < 8}. 

6.2 Particle-hole duality 

A further observation that simplifies the calculation of the eigenfunctions is 
a particle- hole duality result. We introduce a dual operator, in the spirit of 
the Hodge star operator from differential geometry (see e.g. [Jos02j ). by 

Definition The dual of a\^), with 

being any element of the basis / T^Oj) and a G C, denoted by >k(q;|\1')), is given 
by 

*(a|vl>)):=a*a(^,^_J...a(V'n)|l), (34) 

where 

1 1) := IV^ls T V^lsi V52s T V52s i V2p^ T V52pi i V2p2 T 'P2p2 i V2p3 T V52p3 i) 

and a{ip) is the usual annihilation operator which maps \il) ipi^ . . . ipi^) to 

We extend * linearly to real linear combinations of the a|\&), thereby obtain- 
ing an antilinear map from Vo(A^) to Vo(10 — (A^ — 2)). We then have the 
following result: 
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Lemma 6.2. Suppose "^ e Vo(N) satisfies L^"^ = £\1/ and S^\l/ = S"^ . 

Then L^(*^) = /:(*^) anc? S^(*^) = 5(*^). Furthermore, if L^^ = M^, 
5*3^ = M^^ and M = p^, then L3(*^) = -M(*^), S'3(*^) = -M^(*^) 
anc? ^(*\l/) =p(*\l/). 

Proof Direct calculations using the second quantized forms of L and S show 
that both operators anticommute with * on Vq{N). The results for angular 
momentum and spin are then trivial. The result for the inversion operator 
follows from the fact that the parity of a wavefunction is equivalent to the 
parity of the number of p-orbitals present in each Slater determinant (since 
R(PnsC( = (PnsO: for both n = 1 and 2, and Rip2p^a = —(p2piOi for z = 1, 2, 3) 
and the number of p-orbitals in the dual of a Slater determinant with k p- 
orbitals is 6 — fc, preserving the parity. D 

The above result is a modest generalization of insights by spectroscopists 
(who did not know they were speaking what mathematicians would call su- 
persymmetry) . They introduced the dual of a configuration with respect to a 
single open shell (termed "conjugate configuration" in |Con80] ). and noticed 
that it gives rise to the same set of L3 and S3 eigenvalues [ConSOj and, what 
is more, the same L^, S^, and L ■ S matrices |CS35j . 

6.3 Simultaneous L -S Eigenspaces 

We now form the joint angular momentum and spin eigenspaces within 
Vo(iV). Lemma [6.21 shows that we only need to do this for Lithium-Carbon, 
the remaining cases follow by duality. Under the restriction of a fixed num- 
ber of Is, 2s, 2p orbitals, the results were no doubt known to early spec- 
troscopists, who realized that the multiplet structure of observed spectra 
can only be captured via superposition of aufbau principle Slater deter- 
minants into "terms" (in our language, joint L^-S^ eigenspaces); see e.g. 
|CS35| IConSO] . We do not however know of a complete tabulation. 

Theorem 6.1. For the Lithium-Neon sequences (N = 3, . . . , 10, Z > 0), 

orthonormal basis for the L^-S^ simultaneous eigenspaces within Vo(A^) are 
as given in Tahles^ni\ See below for the notation used in the tables. 



Proof We only give the proof for the highest dimensional case. Carbon, the 
other cases being analogous but easier. 

By Lemma 16.11 it suffices to find the joint L^-S^-eigenstates in the four- 
electron vector space spanned hj B = {\4'i-^4'i24'i34'ii) \i < ii < i2 < is < 
H < 8}, with the ipi as in ( l29l) . 

We note first that each Slater determinant in the above basis is already 
an eigenfunction of S3, and that the space with S'3-eigenvalue M is isomor- 
phic (by fiipping all spins) to that with eigenvalue —M. Since both L^ and 
S^ commute with 5*3, it suffices therefore to consider their action on the 
eigenspaces of S3, with eigenvalue M > 0. 
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Next we observe that within each such S'3-eigenspace, the span of those 
Slater determinants which share the same number of different spatial orbitals 
(4, 3, or 2) is also invariant under L^ and S^. 

We now calculate the matrices of L^ and S^ with respect to the so- 
obtained subsets of the basis B, using ( l32l) . ( l33l) . and the formulae from 
subsection 16. 1[ To shorten the notation, we will write s t, s |, p, j, Pi I 
instead of (/^s^ T, V^24, f2p, T, f2p, i- 
Four different spatial orbitals, M=2 On |st pi t P2T PsT); 



L' = 0, S^ 



6. 



Four different spatial orbitals, M=l With respect to the basis {|st pi 

P2T Psi), kt Pit P2i Ps]), |sT pii P2T PsT), ki Pit P2T PsT)}, 



/^ 4 


-2 


-2 \ 




/ 3 1 1 1 \ 


-2 
-2 


4 
-2 


-2 
4 


, s^ = 


13 11 
113 1 


V 





y 




VI 1 1 3^ 



Four different spatial orbitals, M=0 With respect to the basis {\s^ pi f 

P2i Psi), kt Pli P2I Psi), \s] Pli P2i PsT), \si Pli P2T PsT), ki 

Pit P2i Pat), |si Pit P2T Psi)}, 

/ 2 1 1 1 1 \ 
12 110 1 
112 110 
112 11 
10 112 1 

\ 1 1 1 1 2 / 

Three different spatial oritals, M=l On each Slater determinant |s| sj. 
Pit PjT)) and each Slater determinant |pfcT Pkl Pit Pit); 

L2 = 2, S^ = 2. 
(In total, these span a 6-dimensional subspace.) With respect to each basis 

(bit Pii sT Pit), bfct Pfci sT Pit)}, 



/ 4 


-2 


-2 








0\ 


-2 


4 


-2 











-2 


-2 


4 




















4 


-2 


-2 











-2 


4 


-2 


V 








-2 


-2 


4/ 



4 
-2 



(In total, these span a 6-dimensional subspace.) 

Three different spatial orbitals, M=0 With respect to each of the bases 

{|sT si Pi] Pji), |st si Pii pjt)} and {Ipfct Pk i Pi] Pj i), \Pk] Pk i 

Pii Pit)}, 

' 4 -2 \ ^2 ^ /^ 4 -2 

-2 4 r - 1-2 4 
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(In total, these span a 12-diniensional subspace.) With respect to each of 
the bases {|piT Pii s] Pji), \pk] Pkl s] Pj[), \pi] pd s[ pj]), \pk] 
Pki si Pj\), 



1 4 


-2 





0\ 




/ 1 1 \ 


-2 



4 





4 



-2 


, s^ = 


10 1 
10 10 


V 





-2 


W 




^0101/ 



(In total, these span a 12- dimensional subspace.) 

Two different spatial orbitals, M=0 With respect to the bases {\s] s[ 
Pi T Pi i), \s] s i p2 T P2 i), \s] s i p3 t P3 i)} and {\pi ] pii pa T P2 i 
), IP2T P2i PsT Psi), IpsT Pai Pit Pii)}, 



4 
-2 
-2 



-2 

4 

-2 



0. 



This completes the explicit description of the action of L^ and S^. 

The eigenfunctions and eigenvalues are now found by explicit diagonal- 
ization of the above matrices. D 

Tables HIITT] use the following conventions: 

1) The two Is orbitals present in every Slater determinant are not shown. 

2) The eigenfunctions are not normalized. 

3) In all cases, it is assumed that i = 1,2,3, {i,j) = (1,2), (2,3), (3, 1) and 
{i,j, k) is any cyclic permutation of (1, 2, 3). In particular, any eigenfunction 
containing a variable corresponds to a 3D subspace. 

4) Eigenfunctions of the form a^i + h^2 + c^z are such that a + h + c = and 
stand for two linearly independent orthogonal choices of {a,b,c), and hence 
correspond to a two-dimensional subspace. 

5) Within each L^-S^-eigenspace, the different S3 eigenspaces are separated 
by a line, in the order M = S, —S, S — 1, —{S — 1), . . . , 0. 

6) The spin orbitals ip2s T, </'2s i, V2p, T, ^2p, i are abbreviated s, s, p,, pi. 

Note that the parity of the eigenfunctions in the tables, although not 
shown explicitly, can be read off by counting the number of p orbitals (see 
the previous section). 

Inspecting these eigenspaces reveals a number of interesting properties. 

Corollary 6.1. For any N = 3, . . . , 10, the maximum dim,ension of any 
simultaneous eigenspace within Vo(A^) of the operators ^ is two. 

As regards diagonalization of the Hamiltonian, this is clearly much more 
promising than the 70-dimensional space of Carbon. 
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s^ = ! 


L2 = 


1^) 


L2 = 2 


\p.) 





8^ = 


8^ = 2 


I? = Q 


\piPi} + \P2P2} + \PiP3} 




V = 2 


\spi) - \sPi) 


\sPi) \PiP3) 


\sPi) \PiPj) 


\sp^) + \spi) \piP]) + \p~iPj) 


V = 6 


\piPj) - IWiPj) 
a\PiPi) + b\P2P2) + cIpsPs) 





Table 4: Lithium 
sequence L^-S^ 

^ P ■ Table 5: Beryllium sequence L^-S^ eigenspaces. 





s^ = ! 


s' = f 


L2 = 


IsPiPl) + \sp2lh) + IspsKl) 


blP2P3) 




IP1P2P3; 


\spiPi) + \sp2P2) + Ispspj) 


|PlP2P3> + blP2P3> + IPTP2P3) 
IP1P2P3) + IP1P2P3) + IP1P2P3) 


L2 = 2 


\ssPi) 

\PiP}Pj) + \PiPkPk) 

2\spiPj) - \sp-pi) - \sPiP]) 

\sspi) 

\p~iPjPj) + \p~iPkPk) 
2\spiPj) - \spiPj) - \spipj) 


\3PiPj) 


\sPiPj) 


\spiP]) + \splPj) + WiPj) 
\sPiPj) + \spiPj) + \sPiPj) 


L2 = 6 


\spip-) - \sp-pj) 
ffl|spipT> + 6|sp2Pi) + clsp-sp^) 

\pipjpj) - \pipkpk) 

a.\p3PiP2) + h\p-iPiP2) + c^PijhPi) 

|sPiP7) - IWlPj) 

ffl|sPiPT> + & SP2PJ) + c\-sp-M) 

\plPjPj) - \plPkPk) 

a\p-iPiP2} + h\p-iPiP2} + C|P3P1P2) 





Table 6: Boron sequence L^-S^ eigenspaces. 



Corollary 6.2. For any N = 3, . . . , 10, and any simultaneous L^-S^-La- 
S3-R eigenspace within Vo(A^) with L3 eigenvalue equal to zero, all Slater 
determinants occuring within the space differ by an even number of orbitals. 

This is remarkable, and will greatly simplify the structure of the Hamil- 
tonian matrix in the basis ( l30l) . due to the simpler structure of Slater's rules 
(see below). Also, it implies that even for the correlated eigenstates of the 
Hamiltonian, the orbitals fl29p are natural orbitals in the sense of Lowdin. A 
more abstract proof of Corollary 16.21 will be given elsewhere. 
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S^ = 


S^ = 2 


S' = 6 


L^ = 


\sspiPi) + \ssp2P2) + sspspj) 

biprP2Pj> + bipTPspj) + |p2PiP3p?> 


3|spiP2P3) - |.SPlP2l>i) - \splP2P'i) - l-SPlPJPs) 


\sp1P2P3} 




\sP1P2P3) 




\sp1P2P3) + \sp1P2P3} + \sp1P2P3} + \sp1P2P3} 


•:ikPiP2P3) - \spiP2P:i) - ISP1P2P3) - ISP1P2P3) 




I.5P1P2P3) + I.5P1P2P3) + ISP1P2P3) 

-|^P2P3> - [spiPiP-i) - \'SPlP2pt) 


\sp1P2P3) + \sp1P2P3) + \sp1P2P3) + \sp1P2P3) 


\sp1P2P3) + \sp1P2P3) + \sp1P2P3) 

+ \WIP2P3) + \SPMP3) + \SP1P2P3) 


L^ = 2 


\sF,PiPj} - ImPjPj) + \spipm) - \sp,pm) 


\sspiPj) 

\PlP2P3Pi) 
\sPiPjPj) + \sPzPkPk) 




\sspipj) 

\PlP2P3Pi) 
IWiPjPj) + WlPkPk) 


\sspip-) + \ssplPj) 

\pipipjpk) + \pipipjpk) 

\splPjPj) + WiPjP]) + \splPkW) + \spiPkPk) 


L^ = 6 


\s'sPiPj) - \ssp-pj) 
a\sspip[) + b\ssp2Pi) + c\ssp3jE) 

\pipipjpk) - \pipipjpk) 

a\PlPlP2P2) + b\PlPlP3P3) + c\P2P2P3Pi) 

\sPlP]Pj) - ^SPiPjPj) - IspIPfcPfc) + IspiPfcpF) 

<i{\sPlP2P3) -+ \SP1P2P3)) + 6(|spiP2P3) + \sPlP2P-i)) 

+c(\splp2pi) + jspipipa)) 


\sPiPjPj) - \spiPkPk) 

<A^PlP2P3) + &|SP1P2P3> + C\SPIP2P3) 




\spiPjPj) - \sp,PkPk} 


a\spiP2P3) + b\spiP2P3) + c\spiP2P3) 


UF.PjP]) + \sp,PjP]) - \splPkPkl - \sPiPkP^) 
a(,\spiP2P3) - \sp1P2P3y) + b(\-'>PiP2P3) - \sp1P2P3)) 

+c{\sp{p2pE) - \sp1IhP3)) 



i^ = ! 



\-VhPiP2P2) + IspiPip-iPa) + \sp2P2PaPa) 



\sP1P1P2P2) + \sp1P1P3P3) + \sp2P2P3P3) 



S'='4 



\ssp1p2p3) 







\SSP1P2P3) 




\SSP1P2P3) + \SSPIP2P3) + 


SSP1P2P3) 


\SSPIP2P3) + \SSPIP2P3) + 


^ 


SSP1P2P3) 


\splp2pipt) 


^ 




\SP1P2P3P^) 




\sPiP^PJPk) + \sp^P^PJPk) + 


\sp^P^PjPk} 



\sspiPjp]) + \ssp,pkPk) 

\plp2p3PJPk) 

2\sp,p[pjPk) - \spip^pk} - \spiplPjPki 
\sWiPjPj) + \ssplPkPk) 

\plp2p3PJPk) 

'2'\spiPtPjPk) - \sptPlPjPk) - \spiPiP]Pk) 



\spiPiP]Pk) + \spiPiPjPk) + \spiPiPjPk) 



\sspiPjPj) - \sspiPkP^) 

a\sspiP2pt) + b\sspipip3) + c\ssplP2P3) 

\sPtP'iPjPk) - \sPtKPjPk) 

a\spiThP2P2) + b\spiplp3pi) + c\sp2THP3P3) 

\ssplP3P3) - \sWiPkPk) 
a\ssp^Pz) + 6|sspipipi) + c\ssplp2pt) 

IspiPiPjPk) - \sPtPiPjPk) 
a\sPlPTP2P2) + b\spiPlP3P3) + C\SP2P^P3P3) 



Table 8: Nitrogen sequence L -S eigenspaces. 





8^ = 


8^ = 2 


L' 


= 


\sSpip[p2P^) + [sspiplpsp^) + \sSp2Pi'P3P3) 
\pmP2mP3P3) 




L' 


= 2 




\ssp1p2P3Pi) 
\sp1P2P3PjPk} 




\ssp1p2p3Pi) 


\spiP2P3PjPk) - \spiP2P3PjPk) 


\sp1P2P3PjPk) 


\sSpiPiPjPk) + \sSPiPiPjPk) 

\splp2p3pjpk} + \splp2p3pjpk) 


L' 


= 6 


\sSPiPiPjPk) - \sSPiPiPjPj) 

a\sspiplp2n) + b\ssp{plp3P^) + c\ssp2piP3P3) 





Table 9: Oxygen sequence L -S eigenspaces 
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s^ = ! 


L2 = 


\sp1PlP2P2P3P3) 


\sp1PlP2P2P3P3) 


L2 = 2 


\ssp1P2P3PiPj) 




\ssp1p2P3PiPj) 





8^ = 


L2 = 


S'spiPlP2P2P3Pt) 



Table 11: Neon sequence L^-S^ 
eigenspaces. 



Table 10: Fluorine sequence L^-S^ 
eigenspaces. 

6.4 Symbolic interaction matrix 

In order to calculate the Hamiltonian matrix PHP on each L^-S^-i? eigenspace, 
note first that PHP commutes with ([7]) (see Lemma HTTj which remains valied 
for the more general orbitals flHTl) . cf. the calculus in Section [HH]). Hence it 
suffices to pick arbitrary components of L and S, say L3 and 5*3, and calculate 
this matrix on the l^-L^-^-S^-R eigenspace with maximal 5*3 and L3 = 0. 
These spaces are shown in Tables [T2HT41 Here, because of their importance 
for the interaction energy, the Is orbitals are shown and the eigenfunctions 
are normalized. We find it convenient to abbreviate the spin orbitals 

V^lsT, ^Isi, ^2s], ^2si, ^2poJ, V52p3i, V52piT, ^2pA^ ^2p2], ^2p2 i (35) 

(even more drastically than in Tables HHH]) by 

1, T, 2, 2, 3, 3, 4, 4, 5, 5. (36) 

Thus, for example, the top Carbon state of Table [131 

1 



v^ 



(1112233) + 1112244) + |112255)), 



stands for 

\Vls^ Vlsi V2sl V2si V2P3I '^2pd) + IV^lsT Vlsi V2s^ <^2si <^2pil V2pii) 

+ Iv^lsT flsi f2s] f2si f2p2l f2p2i)]- 



73 V 



We begin by analyzing the Vee matrix elements between the eigenfunctions 
of Tables [T2l - [T4l Using Slater's rules |SU96[ Section 2.3], these are straight- 
forward to express in terms of Coulomb and exchange integrals {aa\bb) and 
{ab\ba) of the spatial orbitals (IHTj) . where - in common notation - 



{ab\cd) 



dxidx2a* {xi)b{xi) - 



1 

\Xi — X2I 



rC*{x2)d{x2). 



(37) 
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Lemma 6.3. Let A^ G {3, . . . , 10}. Orthonormal bases of the simultaneous 
L^-^-L^-S-i-R eigenspaces within Vo{N) with S3, maximal and L3 = 0, and 
the corresponding V^e matrix elements (\l/|Vee|^E') in terms of Coulomb and 
exchange integrals of the one- electron orbitals [3l\}. are as given in Tables 
[TB^4\ Here the orbitals are abbreviated as in [3^ } - [3^) . and the off-diagonal 
matrix elements ("$ 7^ \E'j in the two-dimensional eigenspaces are denoted by 
"cross". 

The shortness of the expressions for the (\l/|\4e|^); and the absence of 
Coulomb and exchange integrals involving the last orbital, comes from the 
absence of single excitations (Cor. 16.21) and the equivalence of the p orbitals 
in ( pTl) up to rotation. The latter would be destroyed by changing to a 
basis of L3 eigenfunctions (which is why we have not done so even though 
this would have been more convenient for the diagonalization of L^ in the 
previous subsection). 

6.5 Explicit interaction matrix 

In order to obtain explicit values, we finally need to substitute the explicit 
PT orbitals ( 12^ . ( I25ll . and evaluate the ensuing Coulomb and exchange in- 
tegrals. We do this via a four-step procedure: reduce the original integrals 
over R^ to integrals over M^ via Fourier transform calculus; explicitly de- 
termine the Fourier transforms of pointwise products of the above orbitals; 
reduce to ID integrals with the help of spherical polar coordinates in Fourier 
space; evaluate the remaining ID integrals, which turn out to have rational 
integrands. 

The Fourier transform of a function / G L^(R") will be denoted /; we 
find it convenient to use the definition 

fik) := / f{x)e-'^-^dx (38) 

which does not contain any normalization constants. 

Lemma 6.4. For one-electron orbitals ipa G L^(]R^) fl L°°(]R^), with ip^ ^ 
L2(M3) n L°°(R3) and a G {i,i,k,l}, let f{x) := ^i{x)ip*{x) and g{x) := 
ipl{x)ilJe{x) . Then 

ii^iipjlipki^e) = / dxidx2tpi{xi)tljj{xi)- 7ipl{x2)i^e{x2) 

dk-^Arr{k)g{k). (39) 



27r2 7j,3 |A;|2 

Note that this shows that exchange integrals [ipiipjlipjipi) are positive. 

Proof Since the tpa are in L^(M^) fl L°°(M^), their products / and g are 
in L^(M^) n L°°(R^) and hence their Fourier transforms are well defined. 
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1/ 


s' 


i? 


* 


(Ke) 


Li 





3 

■1 


1 


|1T2) 


(11|11) + 2(11|22)-(12|21) 


2 


3 

■1 


-1 


1113) 


(11|11) + 2(11|33)-(13|31) 


Be 








1 


1T22) 


(11|11) + 4(11|22) - 2(12|21) + (22|22) 


^(|1T33) + |1T44) + |1T55)) 


(11111) + 4(11|33) - 2(13|31) + (33|33) + 2(34|43) 


cross 


v/3(23t32) 


2 





-1 


^(|1T23)-|1T23)) 


(11|11) + 2(11|22) - (12|21) + 2(11|33) - (13|31) 
+ (22133) + (23|32) 


2 


2 


-1 


lT23) 


(11111) + 2(11|22) - (12|21) + 2(11|33) - (13|31) 
+ (22|33) - (23|32) 


1 


1T45) 


(11|11) + 4(11133) - 2(13131) + (33144) - (34143) 


6 





1 


^(2|lT33)-|lT44)-|lT55)) 


(11111) + 4(11133) - 2(13131) + (33133) - (34143) 


B 





3 

4 


1 


^ (|1T233) + |lT244) + |lT255)) 


(11|11) + 2(11|22)-(12|21) + 4(11|33)- 2(13131) 
+2(22|33) - (23|32) + (33|33) + 2(34|43) 





15 
1 


-1 


1T345) 


(lljll) + 6(11133) - 3(13131) + 3(33144) - 3(34143) 


2 


3 

4 


-1 


1T223) 


(lllll) + 4(11122) - 2(12121) + 2(11133) - (13131) 
+ (22122) + 2(22133) -(23132) 


^(|1T344) + |1T355)) 


(lllll) + 6(11133) - 3(13131) + (33133) + 2(33144) 


cross 


y2(23|32) 


1 


^ (2111245) -|1T245)-|1T245)) 


(lljll) + 2(11|22) - (12|21) + 4(11|33) - 2(13|31) 
+2(22133) + (23|32) + (33|44) - (34|43) 


2 


ir, 
1 


1 


1T245) 


(lllll) + 2(11|22) - (12|21) + 4(11|33) - 2(13|31) 
+2(22|33) - 2(23|32) + (33|44) - (34|43) 


6 


3 
4 


1 


^ (2|lT233) - |lT244) - |lT255)) 


(lllll) + 2(11122) - (12121) + 4(11|33) - 2(13131) 
+2(22|33) - (23|32) + (33|33) - (34|43) 


-1 


^ (2|lT345) - |1T345) - |lT345)) 


(lljll) + 6(llj33) - 3(13j31) + 3(33j44) 



Table 12: Vee matrix element expressions for the Li-B sequences, 'cross' de- 
notes the off-diagonal term in the 2x2 matrix. See Lemma [6.31 for notation. 
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L^ 


S' 


fl 


* 


(Vee) 


c 








1 


^ (|1T2233) + |1T2244) + |lT2255)) 


(11|11) + 4(11|22) - 2(12|21) + 4(11|33) - 2(13|31) 
+(22|22) + 4(22|33) - 2(23|32) + (33|33) + 2(34|43) 


^ (|1T3344) + |1T3355) + |lT4455)) 


(11|11) + 8(11|33) - 4(13|31) + 2(33|33) + 4(33|44) 


cross 


2(23|32) 





2 


-1 


^(3|1T2345) - I1T2345) 

-1112345) -1112345)) 


(11|11) + 2(11|22) - (12|21) + 6(11|33) - 3(13|31) 
+3(22|33) + (23|32) + 3(33|44) - 3(34|43) 





6 


-1 


|1T2345) 


(11|11) + 2(11|22)-(12|21) + 6(11|33)-3(13|31) 
+3(22|33) - 3(23|32) + 3(33|44) - 3(34|43) 


2 





-1 


i(|lT2344) -1112344) 

+ |1T2355)-|1T2355)) 


(11|11) + 2(11|22) - (12|21) + 6(11|33) - 3(13|31) 
+3(22|33) + (33|33) + 2(33|44) 


2 


2 


1 


1T2245) 


(11|11) + 4(11|22) - 2(12|21) + 4(11|33) - 2(13|31) 
+(22|22) + 4(22|33) - 2(23|32) + (33|44) - (34|43) 


|1T3345) 


(11|11) + 8(11|33) - 4(13|31) + (33|33) + 5(33|44) 
-3(34|43) 


cross 


(23|32) 


-1 


^(|lT2344) + |lT2355)) 


(11|11) + 2(11|22) - (12|21) + 6(11|33) - 3(13|31) 
+3(22|33) - 2(23|32) + (33|33) + 2(33|44) 


6 





1 


^ (2|lT2233) - |1T2244) - |lT2255)) 


(11|11) + 4(11|22) - 2(12|21) + 4(11|33) - 2(13|31) 
+(22|22) + 4(22|33) - 2(23|32) + (33|33) - (34|43) 


^ (2|lT4455) - |1T3344) - |lT3355)) 


(11|11) + 8(11|33) - 4(13|31) + 2(33|33) + 4(33|44) 
-3(34|43) 


cross 


-(23|32) 


-1 


^(2|lT2345) - |1T2345) - |lT2345) 
+2| 112345) - 1112345) - 1 112345)) 


(11|11) + 2(11|22) - (12|21) + 6(11|33) - 3(13|31) 
+3(22|33) + 3(33|44) 


6 


2 


-1 


Jg(2|lT2345) - |1T2345) - |lT2345)) 


(11|11) + 2(11|22) - (12|21) + 6(11|33) - 3(13|31) 
+3(22|33) - 2(23|32) + 3(33|44) 



Table 13: Vee matrix element expressions for the C sequence, 'cross' denotes 
the off-diagonal term in the 2x2 matrix. See Lemma 16.31 for notation. 
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L^ 


s^ 


R 


* 


{Vee) 


N 





3 

4 


1 


^(|lT23344) + |lT23355) 
+ |1T24455» 


(11(11) + 2(11122) - (12121) + 8(11133) - 4(13|31) 
+4(22133) - 2(23132) + 2(33133) + 4(33144) 





13 
1 


-1 


|1T22345> 


(lllll) + 4(11122) - 2(12121) + 6(11133) - 3(13131) 
+ (22(22) + 6(22133) - 3(23132) + 3(33144) - 3(34143) 


2 


3 
1 


-1 


1T34455) 


(11(11) + 10(11(33) - 5(13(31) + 2(33133) + 8(33(44) 
-4(34(43) 


^(|lT22344) + |lT22355)) 


(11(11) + 4(11(22) - 2(12121) + 6(11133) - 3(13131) 
+ (22(22) + 6(22(33) - 3(23(32) + (33(33) + 2(33(44) 


cross 


\/2(23(32) 


1 


^(2|lT23345)-| 1123345) 
-[1T23345)) 


(11(11) + 2(11(22) - (12(21) + 8(11(33) - 4(13(31) 
+4(22(33) + (33(33) + 5(33(44) - 3(34(43) 


2 


ir, 
1 


1 


1T23345) 


(11(11) + 2(11(22) - (12(21) + 8(11(33) - 4(13(31) 
+4(22(33) - 3(23(32) + (33(33) + 5(33(44) - 3(34(43) 


6 


3 

4 


-1 


^(2llT22345)-|lT22345) 
-|1T22345)) 


(11(11) + 4(11(22) - 2(12(21) + 6(11(33) - 3(13(31) 
+(22(22) + 6(22(33) - 3(23(32) + 3(33(44) 


1 


^(2|lT24455)-|lT23344) 
-|1T23355)) 


(11(11) + 2(11(22) - (12(21) + 8(11(33) - 4(13(31) 
+4(22(33) - 2(23(32) + 2(33(33) + 4(33(44) - 3(34(43) 











1 


^(|1T223344) + |1T223355) 
+|1T224455)) 


(11(11) + 4(11(22) - 2(12(21) + 8(11(33) - 4(13(31) 
+ (22(22) + 8(22(33) - 4(23(32) + 2(33(33) + 4(33(44) 


|1T334455) 


(11(11) + 12(11(33) - 6(13(31) + 3(33(33) + 12(33(44) 
-6(34(43) 


cross 


%/3(23(32) 


2 





-1 


^ (|1T234455) - I1T234455)) 


(11(11) + 2(11(22) - (12(21) + 10(11(33) - 5(13(31) 
+5(22(33) - (23(32) + 2(33(33) + 8(33(44) - 4(34(43) 


2 


2 


-1 


|lT234455) 


(11(11) + 2(11(22) - (12(21) + 10(11(33) - 5(13(31) 
+5(22(33) - 3(23(32) + 2(33(33) + 8(33(44) - 4(34(43) 


1 


|lT223345) 


(11(11) + 4(11(22) - 2(12(21) + 8(11(33) - 4(13(31) 
+ (22(22) + 8(22(33) - 4(23(32) + (33(33) + 5(33(44) 
-3(34(43) 


6 





1 


^(2|lT224455) - |lT223344) 
-|1T223355)) 


(11(11) + 4(11(22) - 2(12(21) + 8(11(33) - 4(13(31) 
+ (22(22) + 8(22(33) - 4(23(32) + 2(33(33) + 4(33(44) 
-3(34(43) 


F 





3 
4 


1 


|1T2334455) 


(11(11) + 2(11(22) - (12(21) + 12(11(33) - 6(13(31) 
+6(22(33) - 3(23(32) + 3(33(33) + 12(33(44) - 6(34(43) 


2 


3 
4 


-1 


1T2234455) 


(11(11) + 4(11(22) - 2(12(21) + 10(11(33) - 5(13(31) 
+(22(22) + 10(22(33) - 5(23(32) + 2(33(33) + 8(33(44) 
-4(34(43) 


Ne 








1 


1T22334455) 


(11(11) + 4(11(22) - 2(12(21) + 12(11(33) - 6(13(31) 
+ (22(22) + 12(22(33) - 6(23(32) + 3(33(33) + 12(33(44) 
-6(34(43) 



Table 14: Vee matrix element expressions for the N-Ne sequences, 'cross' 
denotes the off-diagonal term in the 2x2 matrix. See Lemma [U^ for notation. 
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Considering the integral /(A) = ^3 J^3 dxdy " |^^^|" f*{x)g{y), A > 0. It is 



easy to show that 



47r 



Since f,g e L\W) n L 



'oo/]n)3\ 



\^-v\ 



\x\ A2+|fc| 

f*{x)g{y)\ G L^(]R^) and so by dominated convergence 



it follows 



that i^^^^^^*''--'!-'''"'!! c rWraje 



Setting h = e"^!^!; 



dxdy- 



1 



F-yl 



-r(x)^(y) (A^o). 



(40) 



|a;| we have /(A) = j^^ dy{f * h){y)g{y)dy and f * h = 
fh e L1(R3) since / G ^^(R^) and /i G L°°(R3). By Plancherel's theo- 
rem we have /(A) = t^;^ Jj^g (ifc ^a^T^p {f)*{k)g{k), and again by dominated 
convergence and using f,g E L^(R^) fl L°°(R^), 



/(A) 



27r2 



1 

1^' 



^^m^(/)*(^M^) (^^0)- 



Combining ( l40l) and ( 14T1) gives the result. 



(41) 
D 



Next we calculate the Fourier transforms of pointwise products of the hy- 
drogen orbitals ([211), ^^M- Here and below, by expressions such as |x|e~'^l^l(A;) 
we mean the Fourier transform f(k) of the function f{x) = |a;|e~'^I^L 

Lemma 6.5. With the Fourier transform as defined in [3^) . and A > 0, 



Function 


Fourier transform 


g-A|x| 

\x e"^l^l 
\x ^e-'^l^l 

^2g-A|x| 


SAtt 


(A2 + |fc|2)2 
32A2 7r Sir 


(A2 + |fc|2)3 (A2 + |fc|2)2 
192A37r 96A7r 


(A2 + |fc|2)4 (A2 + |fc|2)3 
32i\nkj 


(A2+|fc|2)3 
32A^ 192A7rfc2 


(A2 + |fc|2)3 (A2 + |fc|2)4 
192A7rfej ki 
(A2+|fc|2)4 

32inkj 192i\'^TTkj 


(A2 + |fc|2)3 (A2 + |fc|2)4 



Proof Let f{x) := e~^^^^. We have f{k) = 87rA/(A^ + |A;p)^, which is easy 
to verify by direct calculation (convert to spherical polar coordinates and 
integrate). All remaining Fourier transforms can be expressed in terms of 

_di„-A|x| 



derivatives of/, as follows. Using \x\e "'"' = —^e 



x\x\ ^p~M^\ 



|^|2g-A|x| 



dA2" 



and noting that differentiation with respect to A commutes with the Fourier 
transform gives 



\^^-m = -^m, 



-M-\{k) 



dX- 



■f{k). 
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For the next three Fourier transforms, we recall the well known differentiation 
identities for Fourier transforms: 



Xjf{k) = i-QT.fi^) ^^^ XiXjf{k) 
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dkpdk. 



■fik). 



Consequently 



d 



{xje-^\-\){k)=i—-f{k), xjXee-^\-\{k): 
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dkjdk, 



■m. 



The final case needed is 



{\x\xje-^\''\){k) 



d /, d 
dX \ dk, 



fik) 



Working out the above derivatives of / explicitly is straightforward, yielding 
the formulae given in the lemma. D 

Lemma 6.6. The Fourier transforms of pointwise products of the hydrogen 
orhitals [24^ , ( f^) are as given in the following table. In all cases j, i = 1,2, 3, 
3^t 



Function 


Fourier Transform 


01s01s 
02s 02s 
01s02s 

02p,02p, 
01s02p, 
02s 02p, 
02p,02p, 


l&Z^ 


(4Z2 + |fc|2)2 

2Z4 7^6 1 6Z8 


(Z2 + |fc|2)2 (Z2 + |fc|2)3 1 (Z2 + |fc|2)4 

iV2Z^ 9V2Z^ 


((fZ)2 + |fc|2)2 ((|Z)2 + |fe|2)3 

z6 ez^k'^ 


(Z2 + |fc|2)3 (Z2 + |fc|2)4 

6V2iZ^kj 


((|Z)2 + |fe|2)3 

GZ'^ikj SZ^ikj 
(Z2 + |fc|2)4 (Z2 + |fc|2)3 

GkjkeZ^ 

(Z2 + |fc|2)4 



Proof This is simply an application of the results of Lemma [6.51 
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Finally we use these Fourier transforms, along with the reformulation of 
the Coulomb and exchange integrals from Lemma 16.4^ to derive the explicit 
values of these integrals. 



Lemma 6.7. Using the abbreviated notation 1 



-'Is; 



^2s 



, 3 



^2p3, 



^2p 



., 5 



^2p 



J, the Coulomb and exchange integrals (pTy occuring in 



Tables [TM{14\ with the PT orbitals [24\ ), If2^) are given by 



(lllll) 


(11|22) 


(12|21) 


(22|22) 


(11|33) 


(13|31) 


{22|33) 


(23|32) 


(33|33) 


(33|44) 


(34|43) 


57 


177 


16 7 

729 '^ 


77 7 

512^ 


59 7 

243'^ 


112 7 
6561 ^ 


83 7 
512^ 


512^ 


301 7 
2560^ 


447 7 


27 7 
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Proof We insert the Fourier transforms from Lemma [6.61 into fl39l) . change 
to spherical polar coordinates and integrate. The angular integrals are ele- 
mentary and the final radial integrals, which on account of Lemmas 16.61 and 
16.41 have rational integrands, may be evaluated with Maple; for a truly pen 
and paper method, one can use complex contour integration. D 

Note the insteresting multiscale effect that the exchange integrals are 
much smaller then the Coulomb integrals, by about one order of magnitude. 
Nevertheless, as we will see later, the exchange terms play an important role 
in energy level splitting. 

The table in the above lemma together with Tables [JSHH] completes the 
task of evaluating the matrix PVeeP on Vo(A^). 

6.6 The matrix PHP 

The remaining part PHqP of the Hamiltonian PHP is trivial to determine, 
because the space Vo(A^) is an eigenspace of Hq, with eigenvalue given in 
Lemma 15.11 

Ho = PHoP = Z^{-1 - ^y on VoiN). (42) 

By inspection of fH21) . Lemma [UTTt and Tables [T2HI1] we obtain an interesting 
corollary. 

Corollary 6.3. The matrix of the PT Hamiltonian PHP with respect to the 
basis in Tables [TM{l4\ with the orbitals [24\ ), If25\) is a rational matrix. 



7 Atomic energy levels and eigenstates 

The spectral decomposition of the PT Hamiltonians PHP is almost immedi- 
ate from the block form derived in the previous section , the only remaining 
task being the diagonalization of the 2x2 blocks, which may be done ex- 
plicitly: For (orthonormal) wavefunctions \E'i, \E'2 and Ei := (\l/j|i7|^j), the 
eigenvalues are given by 

El + E2 



A± = ^^ ± V(^^)' + mWeemW (43) 

with corresponding normalized eigenstates 



a/1 + cl V ) 



2/1^ 



(44) 
These formulae together with eq. fH^ . Tables [T^HHl and Lemma [6.71 imme- 
diately yield: 
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Theorem 7.1. For N = 3, ..., 10, and Z > 0, the energy levels of the PT 
Model ( fT^) . [F7\ ) are as given in Tables URUU Each eigenspace has the min- 
imal dimension (2L + 1)(25'+ 1) possible for its spin and angular momentum 
quantum numbers L and S (see Lemma \2. ij) . and the up to normalization 
unique corresponding eigenstate with zero L3 and maximal S3 is as shown 
in the Tables. Moreover the levels and eigenstates in the Tables provide the 
leading order asymptotic terms of the true Schrodinger levels as Z ^ 00, in 
the sense described in Theorem \4.1\ 

Note that the ordering of the PT levels is independent of Z, since the spectral 
gaps are linear in Z. 

8 Comparison with experiment and methods 
in the physics and chemistry hterature 

The analytical eigenvalues and eigenstates derived in the isoelectronic limit 
provide a wealth of information on various quantitites of physical and chemi- 
cal interest, and yield a number of insights into the inner working mechanisms 
of the many-electron Schrodinger equation which are not readily available 
from numerical simulations. 

We will discuss, in turn, the obtained L and S values, ground state di- 
mensions, ground state energies, and spectral orderings. 

8.1 L and S values and the notion of 'group' in the 
periodic table 

The ground states themselves are not accessible from experiment, but their 
spin and angular momentum quantum numbers are. As already mentioned, 
the theoretical values agree with the experimental values in every case (see 
Table [3]), not just for large Z but all the way down to neutral atoms {Z = N), 
capturing the nontrivial dependence on the number of electrons. 

An important theoretical feature of L and 5* values as compared to the 
more familiar semi-empirical concept of "hydrogen orbital configurations" 
is that regardless of the approximations made to predict them in practice, 
they remain well defined in the full Schrodinger equation. See Section [21 
It would therefore be of value to base quantum mechanical explanations of 
the periodic table on numbers such as these. In this context we note that 
L and S values suffice to explain quantum mechanically a large part of the 
notion of "group" in the periodic table. Only five different {L, S) pairs occur 
mathematically for the first 10 atoms, and experimentally for the first 20. 
Now these correspond precisely to group 1 (alkali metals), the union of groups 
2 and 8 (alkaline earth metals and noble gases), the union of groups 3 and 7 
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fD 
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i^ 


Ci^ 




p 
o 



Li 



Be 



Symm. 

2 po 



^5* 

3 DO 
1 po 

3p 

15 



2 po 

2p 
2po 



3p 

3 po 

350 

1 po 
3p 

15 



*2 



*4 
*3 



C'l'o 



C*, 



7ip(*3 + C>I'4) 
*6 
*7 
*1 

*2 
*8 



1 



(*6 + OI-t) 
;*9 + C*lo) 



C-to 



*4 

>I'l2 
*8 
*11 
*3 
*5 

^(^e + Otr) 
/i^j (-1.9 + 0*10) 
7T=f(*i + c*2) 



E 



972 I 5965 -7 

8 ' 5832 

972 I 573977 

8 ' 52488 






1679616 
1363969 7 
839808 
2826353 7 
1679616 
1449605 7 
839808 
14673197 7 
8398080 
1 



(2813231 - 5\/1509308377)Z 



11 72 
11 72 



1679016(2813231 + 5v/1509308377)Z 
gyj|^( 16493659 - V733174301809). 

2006759 7 
839808 
40981549 7 



\Z 



o ooyouo 
11 72 I 40981549 

8 ^ 16796160 - 
11 72 I 4151299 7 

8 ^ 1679616 
_ 11 72 I 8322281 7 

8 ^ 3359232 
11^2 _|_ 706213 y 



_ 11 72 , ,„.,-^., . 
8 ' 279936^ 
1172 , 14301407 7 

8 ' 5598720 

-f ^^ + 67ii464 (16493659 + \/733174301809)^ 



-\z- 
-\z- 
-%z- 
-\z- 
-\z- 
-%z- 
-%z- 

-\z- 
-%z- 
-%z- 



{ 3806107 
1 1119744 
/ 19148633 
V 5598720 
464555 7 
139968 
{ 966289 _ 
1 279936 
4730843 7 
1399680 
1904147 7 
559872 ^ 
9625711 7 
2799360 
961915 7 
279936 
242119 7 
69984 

/- 3806107 I 
U119744 ' 
/ 19148633 
I 5598720 
I 966289 I 



\/221876564389)Z 
- 33g|232 ^221876564389) Z 



3359232 



;V62733275266)Z 



33^V 221876564389) Z 
331432^221876564389)^ 
gy|gjgv/62733275266)Z 



-75549 (2\/T509308377 - 6981)\/3 



ggijg (2^1509308377 + 6981)\/3 



j(V733174301809 - 809747)^2 



t(V733174301809 + 809747)a/2 



-ggij5(V221876564389- 



460642) 
5iTF(\/221876564389 - 460642) 



-5giTr(\/62733275266- 230321) 



9g^(V221876564389 + 460642) 

9^(-V221876564389 - 460642) 

gj^(V62733275266 + 230321) 



E (num.) 

-7.0566 
-6.8444 



c (num.) 



-13.7629 
-13.5034 
-13.2690 
-13.0955 
-13.0112 
-12.8377 



-22.7374 
-22.4273 
-22.1753 
-22.0171 
-21.9878 
-21.7612 
-21.6030 
-21.4629 



-34.4468 
-34.3202 
-34.0859 
-34.1838 
-33.7203 
-33.5938 
-33.3688 
-33.3828 
-33.2422 
-32.7641 
-32.6376 
-32.3943 



-0.2311 



4.3271 



-0.1671 



5.9851 



-0.1056 
0.1056 

-0.2047 



9.4668 
-9.4668 

4.8853 
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Symm. 


* 


£; 


C 


E (num.) 


c (num.) 


N 


450 

2£,o 
2 DO 
4p 

2p 
2po 


*2 


13 172 1 2437421 ^7 
8 "^ 559872 ^ 
13 72 1 24551357 v 
8 ^ "^ 5598720 ^ 




-49.1503 
-48.9288 
-48.8195 
-48.1630 
-47.8103 
-47.5888 
-47.5478 
-46.5905 


-0.1324 
7.5532 


-f 2" + 671^464(30636167 - v/1144203315841)Z 

13 v2 1 7549145 v 
8^ + 1679616^ 
13 v2 1 76337819 v 
8 "•" 16796160 
13 v2 1 3843463 v 
8 ^ "^ 839808 ^ 
13 72 1 15393535 '7 
8 "^ 3359232 ^ 


-3g3^(Vll44203315841 - 1032821)V2 


*6 
*8 
*1 
*5 

^(*4 + C*3) 


-fZ'^ + g;^^(30636167+ V1144203315841)Z 


;^^(V1144203315841 + 1032821)v/2 





3p 
ID 

3 po 
ipo 

IS 


*5 
*6 

7=f(*l+C*2) 
*4 
*3 

7it?(*i+^^-*^) 


7 ^72 1 4754911 7 
4^^ ^" 839808 ^ 
7 -72 1 47726257 ^7 
4^ "^ 8398080 ^ 




-66.7048 
-66.5360 
-66.3421 
-65.3265 
-64.8578 
-63.4984 


-0.1458 
6.8592 


-iZ'' + ^^^(9884485 - V89111336179)Z 

7 72 1 9897607 y 
4^ + 1679616^ 


-235^(2x789111336179 - 572179)V3 


( 72 1 1 ("noo/i/ior; 1 , /OOl 1 mfil 70"! "7 


25j=y;(2V89111336179 + 572179)V3 




F 


2po 

^S 


*2 
*1 


15 -72 1 11982943 v 
8 "•" 1679616 
15 ^72 1 4108267 v 
8 ^ "•" 559872 ^ 




-87.6660 
-85.8342 




Nc 


IS 


*1 


n 72 1 2455271 7 
^^ + 279936 '^ 




-112.2917 





(group 3 metals and halogens), the union of groups 4 and 6 (Carbon group 
and Oxygen group), and group 5 (Nitrogen group). See the table below. 




Be KWI c 



Mg 



L = 
S = 



L = l 

S='/2 



Si 



N 



O 



L=l L=0 

S = 1 S = 3/2 



Ne 



Ar 



Moreover, taking into account the gradients of L or S" with respect to atomic 
number A^ would separate the group 3 metals from the halogens, and the 
Carbon group from the Oxygen group. Note that L and S gradients are 
mathematically analogous to ionization energies, which are gradients of en- 
ergy with respect to A^. 

8.2 Ground state dimensions 

These dimensions are shown in Table [31 They are interesting as they are 
a measure of the 'flexibility' within the ground state, in that they specify 
the number of degrees of freedom which can be varied without affecting 
the energy of the state. This flexibility appears to be curiously unexplored 
in the literature, perhaps in part due to it not being clearly captured by 
the semi-empirical Bohr-Slater picture, the Hartree-Fock approximation, or 
Kohn-Sham density functional theory. 

On a qualitative level, we expect that an atom with a high-dimensional 
ground state will form a wider range of molecules than an atom with a 
similar number of valence electrons but with a lower dimensional ground 
state. This should be true both in terms of molecular geometry (e.g. linear, 
bent, triangular, tetrahedral) as well as in terms of which atoms it will stably 
bond with. We plan to develop this idea in a more mathematical way in a 
future publication. 

8.3 Ground state energies 

The asymptotic ground state energies, despite being theoretically justified 
only for strongly positive ions (see Section Ej), still capture around 90% of 
the experimental (RJK"'"07] energies of neutral atoms. See the following table. 



Atom 


Li 


Be 


B 


C 


N 





F 


Ne 


EpT 

^exp 

Error 


-7.0566 

-7.4779 

5.6% 


-13.7629 

-14.6684 

6.2% 


-22.7374 
-24.6581 

7.8% 


-34.4468 

-37.8558 
9.0% 


-49.1503 

-54.6117 

10.0% 


-66.7048 

-75.1080 

11.2% 


-87.6660 
-99.8060 

12.2% 


-112.2917 

-129.0500 

13.0% 
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8.4 Spectral orderings and Hund's rule 

The spectral orderings of the asymptotic levels are in spectacular agreement 
with the experimental data }Huh93l lRJK"'"07] , even for neutral atoms. For 
the purpose of these comparisons we consider only the experimental states 
attributed to configurations containing only orbitals with n < 2. The results 
differ only by the interchange of two higher levels in Beryllium (^D and ^P) 
and Carbon (^D° and ^S°). 

A key virtue of our exact eigenstates is that they allow to trace the spectral 
gaps to the size of individual Coulomb and exchange integrals. 

As an example of a 2s-2p spectral gap, consider the "^S ground state 
and ^P first excited state of Lithium. Table [13] shows that the gap is given 
by the difference in interaction of the 2p and 2s orbitals with the Is shell, 
[2(11|33) - (13|31)] - [2(11|22) - (12|21)]. 

As an example of energy level splitting between two states with an equal 
number of Is, 2s and 2p orbitals, consider the '^S° ground state and ^D° 
first excited state of Nitrogen. A look at Table [TH reveals that the energy 
difference consists only of the exchange term — 3(34 1 43), which is present in 
the ground state due to the parallel spins of the three p-orbitals, but absent 
in the excited state. 

In a large majority of cases, the theoretical orderings also agree with Hund's 
rules. In fact, many of Hund's rules are rigorous theorems in first order 
perturbation theory and related models, and rely only on the structure of 
the symbolic matrices in Tables [T2] - [T^ not their numerical values. This will 
be discussed elsewhere. 
Let us also describe a 

Counterexample to Hund's rules. Consider the higher Carbon ls^2s2p'^ states. 
Hund's rules would order their energies, regardless of the choice of nuclear 
charge Z, as 

Esgo < Ei£)o < Ezpo < Eigo < El DO < Eipo. 

For large Z this agrees with the PT and experimental orderings. (That the 
latter two agree with each other follows from Theorem 14.11 ) 

But experimentally, at Z = 20 the ^D° singlet and the ^S° triplet are 
observed to cross, see [RJK+07] and Figure [H (This crossing is beautifully 
confirmed by theoretical calculations based on a simple CI model designed 
with the help of our asymptotic findings here, see |FG09j .) In particular, in 
the neutral atom, Z = N, the experimental ordering is 

Ebgo < Esdo < Eipo < El DO < Ezgo < Eipo. 

This is an important example because it shows that it is of value to investigate 
which of Hund's rules can be justified quantum mechanically and which ones 



37 



can not. In this particular case, a closer look shows that there should be 
no universal ordering, neither one way nor the other. The energy difference 
as read off from Table [T3] consists of a 2s-2p positive exchange term and a 
2p-2p negative exchange term 



E: 



350 



El 



D° 



(24|42)-3(34|43), 



(45) 



and so could have either sign, depending on the orbitals. 

Note that this interesting effect is missed when the states under investiga- 
tion are modelled by their aufbau principle configurations. By Hund's rules, 
these are 1 112345) for the singlet and 1 112345) for the triplet. A simple calcu- 
lation shows that the energy difference is then -E350 — Ei^o = — (24|42) < 0, 
which is very far from the correct difference (HSll and incorrectly predicts a 
universal ordering. 



8.5 Spectral gaps 

The asymptotic energy levels, despite their excellent orderings, do not give 
quantitatively useful spectral gaps for neutral atoms. However, in the regime 
of highly positive ions in which they were theoretically justified (see Theorem 
14.11) they beautifully match the experimental gaps, as shown in Figure [H 
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Figure 1: Splitting of noninteracting Carbon ground state energy by electron 
interaction. Lines: Asymptotic Schrodinger levels (this paper); Circles: ex- 
perimental data [RJK"'"07l IMooTOj . For the highest level at Z = 6 and the 
fourth level a.t Z = 18, we were unable to find experimental data. 

For the Carbon series (A^ = 6, Z = 6, 7, 8, ... ), we show the experimental 
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spectral gaps 

E^{N,Z) E,{N,Z) 

(circles) and the perturbation-theoretic spectral gaps 

Ef^{N, Z) E[^{N, Z) _ Ef e['^ 

Z2 Z2 ~ ~Z' ~Z' 

(lines) against ^, with the energy of the lowest level (which shifts with Z) 
having been subtracted for clarity. By Theorem 14.1^ the match between 
PT and Schrodinger energy levels would become even better as Z increases 
further. But beyond the value Z = 20 shown here, the match between 
Schrodinger levels and experiment slowly starts to deviate, due to the onset 
of relativistic effects, whose study lies beyond the scope of the present paper. 

8.6 Overall conclusion 

The principal conclusion of this paper is that the semi-empirical hydrogen 
orbital configurations of atoms developed by Bohr, Hund and Slater have a 
precise mathematical meaning, as asymptotic limits of the true Schrodinger 
ground states for large nuclear charge. (This holds up to certain small but 
interesting corrections, as described in Section ISTTl ) We hope that the limit 
eigenstates calculated here (see Table [2]) can serve as a theoretical alternative 
to semi-empirical discussions of the periodic table in the literature. 

Another use of our findings, as benchmark data for the design and vali- 
dation of computational methods, is explored in a companion paper |FG09j . 
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